Universal lineshape of the Kondo zero-bias anomaly in a quantum dot 
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Encouraged by the recent real-time renormalization group results we carried out a detailed anal- 
ysis ol the nonequilibrium Kondo conductance observed in an InAs nanowire-based quantum dot 
and found them to be in excellent agreement. We show that in a wide range ol bias the Kondo con- 
ductance zero-bias anomaly is scaled by the Kondo temperature to a universal lineshape predicted 
by the numerical study. The lineshape can be approximated by a phenomenological expression ol a 
single argument eVsd/fcsTk. The knowledge ol an analytical expression for the lineshape provides 
an alternative way for estimation ol the Kondo temperature in a real experiment, with no need for 
time consuming temperature dependence measurements ol the linear conductance. 
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In quantum dots an unpaired spin of the localized elec- 
tron may act as a single magnetic impurity, which inter- 
acts with electrons in metallic leads forming a Kondo- 
correlated state The quantum dots have been rec- 
ognized as a versatile system for probing the many-body 
nature of the Kondo effect @, H| ■ The main advantage 
of quantum dot hosting the Kondo state is the ability to 
tune its main characteristic parameter, the Kondo tem- 
perature (Ik). Thereby, it is possible to study the re- 
action of the many-body system to various external per- 
turbations, such as temperature T, bias V and magnetic 
field. 

One of the remarkable properties of the Kondo effect is 
that the response of the Kondo-enhanced conductance to 
a perturbation is governed by a set of universal laws inde- 
pendent of the physical system in which the Kondo state 
is realized. It has been experimentally verified that in 
the low-energy limit ({T, e|V|/fcB} << ?k, where fcs is 
the Boltzmann constant and e is the elementary charge) 
the Kondo conductance as a function of T and V is 
described by a universal quadratic law 0, [B| expected 
from the Fermi-liquid theory @, 0] • At higher tempera- 
tures (T ~ Tk) the Kondo conductance has been shown 
to scale to a universal dependence of a single parame- 
ter T/Tk found from numerical renormalization group 
(NRG) study @ . These NRG calculations were later ap- 
proximated by the phenomenological expression Q of a 
single parameter T/Tk, which became the main method 
for estimation of the Kondo temperature in real systems. 
Recently, we also attempted to check the universal be- 
havior of the magnetic field dependence in both the low- 
and intermediate-field range [5|. Despite all the theoret- 
ical advances, the problem of the nonequilibrium Kondo 
model at intermediate bias (|V| ~ k^T^/e) remained 
unsolved. As a result, the lineshape of one of the main 
hallmarks of the Kondo effect in quantum dots, known as 
zero-bias conductance anomaly (ZBA) H had not been 
fully described by the theory. Only very recently, the 2- 
loop real-time renomalization group (RTRG) calculations 
developed for the Kondo model [9( provided a numerical 



description of the nonequilibrium Kondo conductance in 
a wide range of biases. The result of these calculations 
agreed very well with the experimental data and linked 
previously developed analytical theories made in the low- 
j(| and high-energy [l(| limits. The development of the 
above mentioned RTRG calculations gave us an addi- 
tional impulse to study the Kondo ZBA in more detail 
and test the universality of its lineshape. 

In this Letter we present a detailed study of ZBA 
associated with the Kondo effect observed in an InAs 
nanowire-based quantum dot. The analysis of our ex- 
perimental data shows that in a wide range of applied 
bias the Kondo ZBA can be scaled into a universal line- 
shape predicted by the RTRG calculations [§] . We found 
our experiment to be well described by the proposed 
phenomenological expression, which accurately approx- 
imates the RTRG results in the experimentally relevant 
range of bias. We also show that this expression can be 
employed as an alternative method of quick and accurate 
estimation of the Kondo temperature. 

Experiment. -The quantum dot used in the experi- 
ment was formed in a 50 nm-diametcr high-quality InAs 
nanowire grown by the vapor-liquid-solid method on a 
(Oil) InAs substrate using molecular-beam epitaxy. The 
as-grown nanowires were randomly deposited onto a p + - 
Si/Si02 substrate and individually connected by ohmic 
(Ni/Au) source and drain electrodes using e-beam lithog- 
raphy. To avoid undesirable effects of the underlying 
substrate the nanowires were suspended in vacuum over 
grooves predefined in the substrate, and fixed from the 
sides by the contacts. The lateral size of the dot was 
defined by the contacts separation (~450 nm) due to the 
electrons in nanowire being localized between two Schot- 
tky barriers formed at the nanowire-contact interface. 
The experiment was performed in a He 3 -Hc 4 dilution 
refrigerator with the base temperature Tb aS c ~10 mK. 
The transport measurements were made by a standard 
lock-in technique. Depending on the temperature the ac 
excitation bias was kept smaller or equal to ksT. The 
dc bias V s d was applied to the source with respect to the 
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FIG. 1. (a) The temperature dependence of the linear con- 
ductance observed in the InAs nanowire-based quantum dot. 
(b) The two-dimensional gray-scale plot of the nonequilibrium 
conductance measured in V g — V 3 d plain at T = Tbase for the 
same range of V g as in (a). The colored vertical dashed lines 
mark the conductance traces at fixed values of V g used for the 
scaling analysis (V g = —2.85 V brown, V g = —2.845 V red, 



-2.84 V orange, V g 



-2.835 V yellow, V g 



-2.83 V 



green, V g 



-2.825 V cyan, V g 



-2.82 V light blue, 



V g = -2.815 V dark blue). 



drain, " virtually" grounded through the transimpedance 
preamplifier. The dot occupancy was tuned by the back- 
gate voltage, V g , applied to the p + -Si substrate. More 
details on growth, fabrication and experimental set-up 
can be found in Ref. [j| 

Large bias scaling.-To investigate the Kondo ZBA in 
more details we used the experimental results previously 
examined in our study of the low-bias and magnetic field 
scaling of the spin- 1/2 Kondo conductance [see Ref. @, 
Fig. 1(c)]. Figure HJa) presents the linear differential con- 
ductance G through the quantum dot as a function of V g 
taken at different T. The region of V g at around —2.83 V 
corresponds to the conductance enhanced at lower tem- 
peratures due to the presence of the many-body Kondo 
state. There are two Coulomb blockade peaks emerg- 
ing at higher temperatures, which mark the region of 
V g with an odd dot occupancy. To identify the relevant 
Kondo ZBA in Fig. QJb) we show the gray-scale plot of 
the nonequilibrium differential conductance measured for 
the same range of V g . The ZBA is seen here as a black 
horizontal line at around V sc i = 0. To test the theo- 
retically predicted scaling of the Kondo ZBA we chose 
several nonequilibrium conductance traces taken at fixed 
values of V g as shown in Fig.QJb) by colored dashed lines. 
For each of the chosen values of V g we estimated 7k by 
fitting the temperature dependence of the linear conduc- 
tance G(T) with the phenomenological expression Q 



G(T)/G = 1 + (T/T£) 



(1) 



where G = G(T = 0,V sd = 0), T K = r K /(2V« - I) 1 / 2 
and the parameter s — 0.22. Here the definition of 7k is 
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FIG. 2. (a) Scaled temperature G(T) dependence measured 
at different V g (colored circles) [see Fig. HJb) for color code], 
approximated with Eq. {TJ (solid black curve), (b) The depen- 
dence of the estimated value of Tk on V g . The error bars here 
represent the 68% confidence interval, (c) The normalized 
conductance G/Go as a function of V s d measure at different 
V g and T — Tb ase . (d) The same normalized conductance 
traces as in (c), but replotted as a function of eVsd/fcsTk 
(colored curves). The traces are shown to collapse onto the 
same universal dependence given by RTRG calculations [s| 
(black dashed curve). 



such that G(T = 7k, V sd = 0) = 1/2G . The experimen- 
tal G(T) dependencies fitted with Eq. ^} with 7k used 
as a fitting parameter are shown in Fig. [Ha). The circles 
represent the experimental points taken at the values of 
V g color coded in Fig. [TJb), and the solid black curve 
corresponds to Eq. ((T|). The deviation from the theoreti- 
cal curve is related to additional mechanisms engaged in 
the transport at higher temperature. To avoid the effect 
of additional mechanisms the fitting procedure was made 
only for T < 200 mK. The fact that at T < 200 mK all the 
experimental data plotted as a function of T/Tk collapse 
on the same theoretical curve reflects the universality of 
the temperature dependence of the Kondo conductance 

||| . The estimated value of Tk at chosen V g is shown 
in Fig. [5][b) and follows the previously observed 0, [ll[ 
parabolic-like dependence [10(. 

As noticed in the early experiments 0, [III Ell the 
width of the Kondo ZBA appears to be proportional to 
Tk; the same qualitative behavior is observed in our ex- 
periment. Figure HJc) shows the normalized nonequi- 
librium conductance G/Gq at different V g with a well- 
pronounced ZBA maximum at V s d = 0. The widest 
ZBA (dark blue curve) is associated with the highest 
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FIG. 3. (a) The RTRG calculations of the nonequilibrium 
Kondo conductance as a function of bias taken from Ref. 
(blue circles) and its approximation with Eq. ((2| (red solid 
curve) and Eq. (O (black dashed curve), (b) Extraction of 
Tk using two different measurements, G(T) (blue squares) 
and G(V) (red circles) fitted with Eq. |T]) (black solid curve) 
and Eq. ([3| (black dashed curve), correspondingly. The two 

(T) (V) 

extracted values and are the same within the sta- 

tistical error, (c) The experimental Kondo ZBA measured at 
V g — -2.83 V [see green dashed line in Fig. []Jb)] (red circles) 
and its theoretical approximation made with Eq. (f3]) (black 
dashed curve). When the bias is such that e|14d| = telk the 
Kondo conductance decreases to 2/3 of its zero- bias value Go- 



Tk ~ 900 mK, and the narrowest ZBA (orange and yel- 
low curves) corresponds to the lowest Tk ~ 300 mK [see 
Fig.[2fb)]. In accordance with the RTRG calculations @ 
the nonequilibrium conductance in the Kondo regime is 
scaled by Tk into a universal dependence (these calcu- 
lations were made in the zero-temperature limit, as for 
our experiment, T/Tk < 0.03, which was theoretically 
checked to be close enough to zero). To verify this 
prediction we plotted the normalized conductance G/Gq 
as a function of the scaled bias eVsd/kBTi in Fig.(2^d). It 
is evident from Fig. [2jd) that all curves collapse onto the 
same universal dependence, which is given by the RTRG 
calculations and plotted as the dashed black curve. The 
deviations from the prediction observed at higher biases 
occur due to the approaching resonant level of the dot. 
In this case the system switches from the Kondo to the 
mixed- valence regime where the scaling is no longer valid 

Phenomenological formula.— As shown above the ZBA 
in the Kondo regime appears to have a universal depen- 
dence on bias, scaled by Tk- This property of ZBA makes 



it potentially valuable for extracting Tk from the bias de- 
pendence of the conductance G(V s d), in the same way as 
for G(T) by approximating it with Eq. (J]). Here the use 
of numerical calculations is quite cumbersome, and the 
authors of Ref. suggested a phenomenological formula 
for describing their RTRG results, which we utilized for 
our definition of Tk 



14] 



G(T = 0,v)/G o 



( 2 V«x - l) v 2 
n + b(\v\" - 1) 



(2) 



where v = eV sd /k B T K . The best fit of Eq. © to the 
RTRG results gives b = 0.05 ± 0.01, si = 0.32 ± 0.01 and 
S2 = 1.28 ± 0.03, Fig. EHa) illustrates the quality of this 
fit. As seen from the plot, Eq. ([2]) (red solid curve) ap- 
proximates the numerical results (blue circles) very well 
in a wide range of bias. However, in real experimental sit- 
uation the Kondo ZBA rarely extends beyond ±2fcBTk/e 
due to a dominating background, thus for the narrower 
range of bias {e\V s d\ / k^T^ < 10) one can use a simplified 
expression 



G(T = 0, v)/Gq = [1 + av 



21 -*i 



(3) 



where a = (2 1 ^ Sl — l)/7r. The simplified formula is iden- 
tical to the one in Ref. [t| As can be seen from Fig. Eta) 
the simplified expression, shown by the black dashed 
curve, approximates the numerical result adequately up 
to eV s d/kBTK ~ 10. Note that in the low-bias limit 
e\V s d\/kBTK -C 1 Eq. ([3]) reduces to the quadratic Fermi- 
liquid dependence G(V s d)/Go w 1 — s\a{eV s d/ k^T^) 2 = 
l — cv(eV s d/kBTK) 2 , where the coefficient cy ~ 0.78 is in 
a good agreement with both theory and experiment 

Since the Kondo ZBA is scaled by Tk into a univer- 
sal lineshape, which is described by a phenomenological 
formula, it should be possible to perform the reverse op- 
eration and extract the value of Tk by fitting the exper- 
imental G(V s d) dependence with Eq. Figure M^>) 
shows a semilogarithmic plot of the Kondo conductance 
measured as a function of T (blue circle) and the modified 
bias eV s d/kB (red circles). Each set of data is approxi- 
mated with its own phenomenological formula using Tk 
as a fitting parameter; G(T) with Eq. Q for T < 200 mK 
and G(Vsd) with Eq. ^ for |V^<j| < 30 /iV, correspond- 
ingly |l5|. The two values of the Kondo temperature 



were found to be T, 



(T) 



314 ± 28 mK for G(T) and 



T K V ; = 293 ± 6 mK for G(V s d) and they are equal within 
the statistical error. Also, fitting of G(V s d) provides sta- 
tistically more accurate value of Tk due to a larger num- 
ber of experimental points available for analysis. In our 
opinion this method of estimating Tk may be advanta- 
geous to the traditional one involving measurement of 
the G(T) dependence and approximating it with Eq. (JlJ. 
First, it is much easier to reliably measure the G(V s d) 
dependence at the lowest possible temperature, rather 
than performing time-consuming measurements of G(V g ) 



4 



in the linear regime at multiple temperatures. Normally, 
the low temperature (T/Tk << 1), required for Eq. (|3]) 
is easily satisfied, especially at submillikelvin tempera- 
tures, since for most systems Tk > 200 mK. Second, as 
pointed out, the value of Tk extracted from the experi- 
mental G(V s d) dependence is potentially more accurate. 

Width of ZBA .-Finally, we would like to discuss an 
important practical implication of Eq. ([3]). So far, the 
width of the Kondo ZBA has been used as a rough ex- 
perimental estimate for Tk @, [HI, [H| and the full width 
at half maximum (FWHM) of the ZBA peak was as- 
sumed to be approximately equal to 2/cbTk- However, 
more careful measurements demonstrated an overesti- 
mation of Tk determined by this method, if compared 
to the values extracted from the G(T) dependence 
To clarify the issue of the ZBA width we did a simple 
mathematical analysis of Eq. ([3]), which revealed that 
G{V sd = ±k B T K /e) = 0.67G « 2/3G and the ZBA 
peak FWHM = 2fc B T K 2y^k B T K Q. By means of 
these simple relations a quick and accurate estimation of 
Tk from the raw data is straightforward. As an exam- 
ple Fig. [3(c) shows the Kondo ZBA peak measured at 
V g = -2.83 V and plotted as a function of the normal- 
ized bias eVgd/fceTK, where Tk is found from fitting of 
the same data with Eq. ([3]) [see Fig. [3(b)]. The universal 
lineshape is shown in Fig. [3(c) by the black dashed curve, 
illustrating an excellent agreement with the experiment. 
It is seen that the width of the Kondo ZBA is 2&bTk 
at two third of its total magnitude, contrary to the ear- 
lier assumed half magnitude [ll|, |l_2j . It also suggests an 
explanation to the discrepancy between the FWHM and 
the value of Tk reported in Ref . HH, alternative to the de- 
phasing by bias [1, [l6| . Unfortunately, we were unable to 
measure the FWHM reliably because in our experiment 
the Kondo ZBA at 0.5Go was already broadened by the 
background. 

In conclusion, using recent results of RTRG calcula- 
tions of the noncquilibrium Kondo model we experimen- 
tally verified the universal scaling of the Kondo conduc- 
tance at intermediate (e|V^| ~ ^bTk) bias. We estab- 
lished that the Kondo ZBA in a quantum dot can be 
scaled by Tk to a universal dependence predicted by 
the numerical calculations and approximated by the phe- 
nomenological formula. An excellent agreement with the 
experiment allowed us to use this formula to extract the 
value of Tk solely from the analysis of the Kondo ZBA 
lineshape at the lowest temperature. This method ap- 
pears to be quicker and more statistically accurate com- 
pared to the traditional one involving the measurement 
of the G(T) dependence. Also, a closer look at the phe- 
nomenological formula revealed that when the applied 
bias is equivalent to fcnTk/e the Kondo conductance is 
at two thirds of its zero-bias value. At the same time 
the FWHM is about 2y / 7rfcBTK, which is larger than it 
was previously thought to be. We demonstrated that 
those relations can provide an immediate and accurate 



estimate of the experimental Tk- 
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